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The Chapman-Enskog method of kinetic theory is applied to two problems of synchronization 
of globally coupled phase oscillators. First, a modified Kuramoto model is obtained in the limit 
of small inertia from a more general model which includes "inertial" effects. Second, a modified 
Chapman-Enskog method is used to derive the amplitude equation for an 0(2) Takens-Bogdanov 
bifurcation corresponding to the tricritical point of the Kuramoto model with a bimodal distribution 
of oscillator natural frequencies. This latter calculation shows that the Chapman-Enskog method is 
a convenient alternative to normal form calculations. 
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PACS numbers: 05.45. +b, 05.20.-y, 05.40.+j 



I. INTRODUCTION 

The Chapman-Enskog method (CEM) has long been 
used to derive hydrodynamic equations of parabolic type 
from kinetic equations of the Boltzmann type pl^| . The 
complexity of these equations has perhaps hindered the 
realization that the CEM is a powerful singular perturba- 
tion method that can be used advantageously as a viable 
alternative to multiple scales or normal form calculations 
in different contexts. Motivated by recent work on the 
subject of models of synchronization of phase oscillators 
|^-^| , 1 shall derive here reduced equations for two differ- 
ent Kuramoto models by means of the CEM. 

The first example corresponds to the limit of small in- 
ertia in a generalized Kuramoto model of globally coupled 
phase oscillators B: 



mujj 



Qj + Kr N sin(^jv - 0j) + 

j = l,...,N, 



(1) 



where 9j, luj and Qj denote phase, frequency and nat- 
ural frequency of the jth oscillator, respectively. The 
natural frequencies are distributed with probability den- 
sity g(£l), which may have a single maximum (unimodal 
distribution), or several peaks (multimodal distribution). 
The positive parameters m and K are the "inertia" and 
the coupling strength, respectively. The complex order 
parameter defined by 
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r N e l ^ N = — V e ie * 
N ^ 



(2) 



measures phase synchronization: in the limit as N — ► oo, 
r N > if the oscillators are synchronized and r N = if 
not. Finally, ^'s are independent identically distributed 
Gaussian white noises, with = 0, (s)) = 

2DSij5(t — s). Without white noise terms, these equa- 
tions were proposed to account for synchronization of 
biological systems jd,||. When the inertial terms vanish, 



m = 0, Eqs. (|l|) and (||) are exactly the usual Kuramoto 
model [p|Jl0[. An important technologically relevant ap- 
plication of this model is the study of superconducting 
Josephson junctions arrays [ pdp^ ]. It has been shown 
that the usual Kuramoto model(with m = 0) describes 
series arrays in the limiting case of weak coupling and 
disorder of arrays with zero junction capacitance pl| . A 
currently open problem is to derive a reduced model for 
the averaged system in the limiting case of weak cou- 
pling, weak disorder and nonzero junction capacitance. 
In so far as capacitance plays the same role as inertia in 
simple electrical networks, the Kuramoto model with in- 
ertia might be useful to understand the effect of nonzero 
electrical capacitance in Josephson junction arrays. 

In the limiting case of infinitely many oscillators, mod- 
els with mean-field coupling are described by an evolu- 
tion equation for the one-oscillator probability density, 
p(6,ijj,n,t), [^3|. For the present model this equation is 



dp 



dp _ D d 2 p 
dt m 2 dup- 

m duj 

where the order parameter is now given by 
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n + jrrBu#-0))p]-«^, (3) 



re 



e lV p(e,uj,n,t)g(n)dndujde. (4) 



Equations (Q) and <Qj should be supplemented with ap- 
propriate initial and boundary data (p is 27r-periodic in 
9 and has suitable decay behavior asw-> ±00) plus the 
normalization condition, 
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+00 



p(9,uj,n,t)dujd9 = 1. 



(5) 



An extense study of oscillator synchronization in the 
model (|) - (§) was carried out in a previous paper Q by 
using different analytical and numerical techniques. The 
limiting case of small inertia, m — ► 0, was analyzed by 
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an arbitrary three-mode truncation of a mode-coupling 
expansion of p 0. It is rather intriguing that exactly 
the same small inertia results had been obtained ear- 
lier by Hong et al Q] using completely different meth- 
ods. Hong et al obtained first a reduced Smoluchowski 
type equation, and then analyzed this reduced equation. 
The derivation of the reduced equation consisted of trun- 
cating a moment hierarchy for the Fokker-Planck equa- 
tion, as explained by Schneider et al |14| ). Both proce- 
dures yielded correctly the linear term of the amplitude 
equation describing the synchronization transition, but 
not the cubic term in this equation. In this paper, we 
shall obtain for the first time a consistent reduced Smolu- 
chowski type equation by using CEM. 

Our second example corresponds to the particular case 
of setting m = and g(fl) = \ [5(tt - O ) + 8 (SI + fi )] 
in (|l|). Then the phase diagram of the incoherent prob- 
ability density p = l/(27r) contains a tricritical point 
where a line of Hopf bifurcations intersects tangentially 
a line of homoclinic orbits at zero frequency. The result- 
ing Takens-Bogdanov bifurcation with 0(2) symmetry 
was analyzed in || by multiple scales techniques, [ p"5[ . 
A shortcoming of these techniques is that all terms in 
the corresponding amplitude equation are necessarily of 
the same order, while terms of different order are needed 
to describe the Takens-Bogdanov point. We shall show 
here how to adapt multiple scale ideas using the CEM to 
overcome this difficulty. Although the resulting ampli- 
tude equation is known, the calculation presented here 
is of independent interest to illustrate the possible ap- 
plications of our presentation of the CEM to deriving 
amplitude and bifurcation equations containing different 
asymptotic orders. 

The rest of the paper is as follows. The Smoluchowski 
equation for the small inertia limit of the Kuramoto 
model with inertia is derived in Section II. Section III 
contains the application of the CEM to finding bifur- 
cation equations. After explaining how the procedure 
works for the well-known Hopf bifurcation, we describe 
how to obtain the normal form corresponding to an 0(2)- 
symmetric Takens-Bogdanov bifurcation for the bimodal 
Kuramoto model (without inertia) . The last Section con- 
tains our conclusions. 



with respect to the corresponding variable). If we call 
ujq a typical unit of velocity, the latter balance yields the 
following velocity scale 




which is just the thermal velocity. The ratio of ojpe to 
m (u>p) u is of the order 



e = V^D, (7) 

which will be our small dimensionless parameter. Lastly, 
we shall choose the time unit so that p t and ujpg are of 
the same order. This choice is dictated by the mechanics 
of CEM (see below) and yields a time unit to = 1/ujq. 
Normalization condition and the definition of the order 
parameter dictate that g(Q) and p are to be measured in 
units of 1/luq too. The angle 9 is already dimensionless. 

Recapitulating, u, fi and K are measured in units of 
the thermal velocity ujq, while p, <?(0) and time are mea- 
sured in units of the reciprocal velocity l/u>o- Then the 
dimensionless Fokker-Planck equation is 

to be solved together with (§), (||), and appropriate pe- 
riodicity, initial and decay conditions as lo — > ±oo. 

A different "parabolic" scaling is usual in kinetic the- 
ory: in it the term [Kr sin(# — ip) — Sl]p u is 0(e) and the 
time derivative is 0(e 2 ). Using the CEM with this scal- 
ing yields the standard Kuramoto model with m = as 
a leading order approximation, while higher order mod- 
ifications contain derivatives of order 4 in 9 and higher. 
We think the "hyperbolic scaling" (ph will yield results 
which are valid over a much larger range of independent 
variables and of parameters. 



A. Chapman-Enskog method 



II. CHAPMAN-ENSKOG METHOD FOR THE 
GENERALIZED KURAMOTO MODEL IN THE 
LIMITING CASE OF SMALL INERTIA 

To find a reduced equation, we should first nondimen- 
sionalize ([}]) and (|J) in an appropriate way. The main 
idea is that the force terms and diffusion in velocity space 
should be dominant. Then the probability density p 
rapidly reaches local equilibrium in velocity and its slowly 
varying amplitude obeys a Smoluchowski equation. This 
means that f2 and K should have the same order as u>, 
and that the terms m~ 1 (ujp) UJ and (D/m 2 )p UJUJ should 
be of the same order (subscripts mean partial derivative 



Setting e = in @, we find a simple equation to be 
solved together with (^) and (||). Its solution is a dis- 
placed Maxwellian: 



V = cj + Kr sin(0 - r/0 - 0, 



p(9,n,t) d9 g (n)dn, 



/ p(9,n,t)d9 = i. 

Jo 



(9) 
(10) 

(11) 
(12) 
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Notice that P(0, ft, t) is an arbitrary function of 9 and t 
except for (|l2|). Furthermore, (||) correspond to a par- 
ticular form of the initial conditions. The Chapman- 
Enskog ansatz consists of assuming that p has the fol- 
lowing asymptotic expansion 



P 



n=l 



pW(6,u>,n;P). (13) 



Furthermore, we impose that the amplitude P obeys an 
equation: 

BP 00 

— = J2e n F^(P), (14) 

where are functionals of P to be determined as the 
procedure goes on. This equation for P is not explicitly 
written in the usual presentations of CEM jl],^. Instead, 
the form of this equation is guessed by writing equations 
for the moments of p and using gradient expansions. We 
find this latter procedure more confusing. 

Insertion of (|13| ) and (|Ti| ) into the equations and aux- 
iliary conditions yields a hierarchy of linear equations 
for the p^ n > . Notice that the latter depend on time only 
through their dependence on P. The functionals F^ n \P) 
are determined so that each equation (and set of auxil- 
iary conditions) for p( n ' has a solution which is bounded 
for all values of u>, even as u) — * ±00. Once a sufficient 
number of F' n ' is determined, ( Jl4] ) is the sought am- 
plitude equation. Please notice that, unlike results from 
the method of multiple scales, terms in ( [Til) may be of 
different order. 

Let us illustrate how the procedure works by finding 
F (0) and Insertion of |jj|) and (|J) in (§), @ and 

(§) yields the following hierarchy of linear equations: 



Cp 



(i) 







Pe 








V 2 
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F (0) 


V2n 







Cp 



(2) _ 



r (0) ^cos(6»- t/>)] 



V2tt 

Ojp { g ] 



o {1) \ 



(15) 



(16) 



and so on. These equations are to be supplemented by 
the normalization conditions ([121) and 



(17) 







p {n) d6duj = 0, 



and the definitions: 

p2ir />oo />oo 

r (n) = 



J — 00 J — 00 



e l(e ~^p {n) d6dLjg{n)dVL, (18) 



Ve~ — 



KA n) sin(9-ip), (19) 



for n = 1, 2, F is now given by (lOl) with r 



,(0) 

given by (11). Lastly, r^ ) and tp are calculated by tak- 
ing the time derivative of ( |Tl| ) and using the first term in 
@ to replace P t by . 

Let us now consider (|l5|). Since the u> integral of its 
left hand side is zero, this equation has a solution only if 
the uj integral of its right hand side is zero. The corre- 
sponding integrals are simplified by using the symmetry 
of the Maxwellian and shifting integration variables from 
w to y. The vanishing of the integral of the right side 
yields 



ir(o) 



{[Xr (0) sin(0 - V) - ft] P] 



(20) 



Notice that we needed F^ ' in the right side of (15) for 
this equation to have an appropriate solution. In turn, 
appeared in this equation due to our choice of the 
time unit. Thus the time unit is dictated by the solvabil- 
ity conditions of the hierarchy of equations generated by 
the CEM. 

Inserting now ( pp| ) in Eq. ([l5]) , we realize that the right 
hand side thereof is the partial derivative of some expres- 
sion with respect to to. This greatly helps finding the 
solution: 



e 2 



2tt 



V 2 



1 



PKr {0) cos(6> - ip) 



+ V 



KP (7 (0) sin(0 - V) - r (o) ^cos(0 - ip) 



- [Fr (0) sin(0 - ip) - ft] r (0) cos(0 - \ - P \, (21) 



which satisfies ( |l7| ) and yields = 0. Notice that a 
term of the form (0) (satisfying J Q w P dQ — 0) could 
have been added to the solution (^l|). However, all 
such terms are already contained in the ansatz ( |l3| ) with 
P = P(6, ft, t; e), and we shall therefore omit them. 

To find F^\ we insert ( pi] ) in ( |l6|) and use the solv- 
ability condition for this equation. Simplifications arise 
from the identities 

/OC O />OC 

-00 J — 00 



and 



(i), d [°° (1) 

O0 J-oo 



The result is 

(i) _ d 



F 



x 



[Pe + ([ifr (0) sin(0 - V) - ft + ij>] 
Kr {0) cos(9 - V) - ^ (0) sin(0 - ^)) f} • 



(22) 



3 



In this expression, f(°) and ip can be found from their 
definition and Eq. (|o|). Integration by parts yields: 

= ^fL _ ( e 3i(fl-^))) +f (n e *(«-*)) 

= tfrW<sin 2 (0-^)) - (fisin(0-^)), 
^ r (o) 

(0 cos(0 - V0> - (sin[2(0 - V)]> , (23) 



r (o) 
where 



/(0,ft)P(0,ft,i;e)d0#(ft)dft. (24) 



We can now insert @, (g|) and @ into (|l|) to obtain 
the sought Smoluchowski equation for P: 

/_K>(0) 
6 — 



Pi 



.r' ' sin( 



- V) - ft] [1 + e^r (0) cos(6> - ip)] 



x JTcos(( 
x Kr {0) (sin : 
+ eP e } = 0. 



sin[2(6» - V)]) - (ftcos(6> 
ip) - eK sm(6 - ip) 



(6-ip)) - (ftsin(0-^>)) 



P 



(25) 



Restoring the original dimensional units, this equation 
becomes: 



P 



d_ 

mK cos 



{{[Kr sin(0 - ip) - ft] [1 + mKr cos(0 - V)] 



( ) (^< s in[2(0 -</,)]> 



- (Clcos(d-ip)))- mK sin(0 - ip) 
x [Kr(sm 2 (8 - ip)) - (fl sin(0 - V))] 
+ £P e } = 0. 



P 



(26) 



to be solved together with (|ll[), definition (|24|), normal- 
ization condition (p"2|), 27r-periodicity in and appropri- 
ate initial conditions. 

Equation (|2^) is the main result of this Section. If 
to = 0, the usual Kuramoto model for phase oscillators is 
recovered. Comparing this result with Hong et al's Equa- 
tion (5), we see that the second and third terms propor- 
tional to to in the drift current were missing and that 
their equation contained an additional diffusive term. 
The latter would be order e 2 in our scheme. These dif- 
ferences can be tracked to the fact that Hong et al fol- 
lowed a previous paper |l4[ ], whose authors used a mo- 
ment method with arbitrary closure assumptions to de- 
rive the Smoluchowski equation. 



III. DERIVATION OF BIFURCATION 
EQUATIONS BY THE CHAPMAN-ENSKOG 
METHOD 

First of all, we shall illustrate the application of the 
CEM to a well-known case: deriving the amplitude equa- 
tion for the usual Hopf bifurcation. Next the much more 



complicated calculation of the amplitude equation for the 
tricritical point of the bimodal Kuramoto model will be 
tackled. 



A. The usual Hopf bifurcation 



Let us consider a system of n ordinary differential equa- 
tions for a n-component vector unknown u{t; a): 



du 

^ = /(u;a); 



(27) 



a is the bifurcation parameter, we shall assume that 
u = is a linearly stable solution if a < and that it is 
unstable if a > 0. All eigenvalues of L(a) = df(0, a)/du 
except A(a) and its complex conjugate have negative real 
parts in a neighborhood of a = 0. We have \(a) ~ 
iujq + ctA'(O) as a — > 0, with RcA'(O) > 0. The linearized 
equation of ( p7j ) about u — 0, a = has the solution 
u = A (p e^M-cc, plus terms which decay exponentially 
fast to zero as time increases. cf>o is the eigenvector cor- 
responding to iu)Q and A is a complex constant. 
To apply CEM to this problem, we shall assume 

a = a.\E + a-2S 2 + 0(e 3 ), 

oo 

u = eA(t; e) <p e^°* + cc + ^ e n u {n) (t; A, A), (28) 

n=2 

where the amplitude A obeys the ansatz 



dA 
~dt 



Y^e n F { - n \A,A). 



(29) 



n=l 



Here A is the complex conjugate of A and cc means com- 
plex conjugate of the preceding term. The parameter e 
measures the size of the bifurcating amplitude and its 
relation to the bifurcation parameter a is determined by 
finding the coefficients ctj. Notice that A varies slowly 
with time (on a scale e 2 t, as we will see). The functions 
vS n > depend on a fast scale t corresponding to stable ex- 
ponentially decaying modes, and on a slow time scale 



through their dependence on A. All terms in (28) which 
decrease exponentially in time will be omitted. a n , 

F (n) 

are determined so that the solutions are bounded as 
t — > oo (on the fast time scale), for fixed A. Substitution 
of (|2S|) and ( p9|) in (p7|) yields the following hierarchy of 
linear equations: 



duV 
dt 



[L(0) - two] 0o = 0, 



- L(0)u< 2 > = [ ax AL'{Q) - P (1) ] fae lu)at + 



cc 



+ -fuu(0:0):(A<p e^ ot + cc) 2 , (30) 
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dt 



L(0> (3) = . .) + [a 2 AL'(0) - F^} 



x 9o e iuiat + cc + -/„„(0; 0) : u {2 \Afoe tbJat + cc) 

+ J/«n M (0; 0) : {Act> e iuot + cc) 3 , (31) 

and so on. Here / uu (0;0) is a n x n matrix, colon sign 
means tensor contraction, etc. The first equation holds 
identically due to the definitions of </>o and iloq. The other 
equations should be solved for bounded as t — > oo. 
Their solutions should not contain terms of the form 
+cc, solving the corresponding linear homoge- 



neous problem. The reason for this latter requirement is 
that all such terms are already contained in A(t;e). 
Eq. © yields 



Oil 



0, F (1) = 0, 



,(2) 



A 2 



[2iwbJ-JK0)]~V«u(0;0):$ 



+cc - |A| 2 L(0)-7««(0; 0) : <^o . (32) 



(/ is the identity matrix of order n) . Solvability of 
yields 



J?W = a 2 X'(0)A- liA\A\ 2 



(33) 



where 



%,9o) 1 



m - ( 3 4, 

(00' 9o) 
2 (^0)/««tt(0;0) : vlfo) 



+<0S,/ UU (O;O) : ^(Oj-VuuCOjO) : <Mo) 



Substituting these results in (|29j), the following ampli- 
tude equation is found: 



^=e 2 a 2 X'(0)A~e 2 fiA\A\ 2 + O(e 3 ). 
dt 



(36) 



These formulas for the bifurcation equation and its co- 
efficients have been obtained many times before; see ]T^ ] 
for an equivalent explicit determination of fi, whose real 
part decides whether the bifurcation is sub or supercrit- 
ical. If Re/i = 0, we should calculate higher order terms 
in the Chapman-Enskog expansion ( |29| ) . The systematic 
way in which the CEM method yields such terms is a 
great advantage with respect to other methods such as 
multiple scales [Ell. 



Similar ideas can be used to derive amplitude equations 
for pattern forming systems governed by partial differ- 
ential equations E^j. In such cases, we have to rescale 
appropriately space variables and assume that the F^ 
also depend on spatial derivatives of A. 



B. Symmetric Takens-Bogdanov bifurcation 

Our starting point in this Subsection is the standard 
Kuramoto model with bimodal natural frequency distri- 
bution, i.e., @ with m = and g(Q) = |[<S(fi - fi ) + 
5(Q + flo)). The phase diagram of the incoherent solu- 
tion P = l/(27r) was depicted in Figure 1 of ||. At the 
tricritical point, P = (K/D = 4,Hq/D = 1) a branch 
of Hopf bifurcations coalesce with a branch with sta- 
tionary bifurcations and a branch of homoclinic orbits, 
in a 0(2)-symmetric Takens-Bogdanov bifurcation point. 
A method of multiple scales was used in || to analyze 
this complicated bifurcation. This method was not com- 
pletely satisfactory because it led to a couple of ampli- 
tude equations whose solutions were later interpreted as 
the two terms of a normal form expansion ||. Let us 
show how a modified CEM leads directly to the ampli- 
tude equation. Near the tricritical point, we may define 
an expansion parameter e so that 



K = K c + K 2 e 2 + 0{e 3 ), 



n = n 0c + n 2 e 2 + o(e 3 ) 

{K C = 4D, n 0c = D). (37) 



The basic slow time scale near the tricritical point is 
T = et, and the method of multiple scales show that 
resonant terms appear in the equations of order e 3 and 
higher Borrowing from the results of that reference 
(which includes making an exponential ansatz for P) , we 
shall make the following Chapman-Enskog ansatz: 



+ 2 m, /««(0; 0) : MHO) - i2Lu I]- L f uu (0; 0) : j . (35) 

Here (j>Q is the eigenvector corresponding to the adjoint 
problem 



P(6,n,t;e) 



1 

2^ 



exp 



£ ^(T i£ ) 

D + in 



J=2 



(38) 



A TT = F<°> (A) + eF^ {A) + 0{e 2 ). (39) 

Terms which decay exponentially on the fast time scale, 
t, will be systematically omitted. The equation for A is 
second order [not first order as (p9|)] because resonant 
terms appear at 0(e 3 ) for the first time. Inserting these 
equations in ( p6| ) (with to = 0), we obtain the following 
hierarchy of linear equations: 

C (a 2 + y) = -±Dd e [axJm e -*<e*\ *i>} 



A T e ie 
D + in 



cc, 



(40) 



02 



d6 = 0, 
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C (a 3 + a x a 2 + y j = -4Dd e {ailm e' ld {e l6 ' , 
<7 2 + f)+(a 2 + f)lme- ie (e i9 ' ) a 1 ) 
+Q 2 Im e-*V 9 Vi)'} 

/ 2 

-# 2 3 e Im e- w <e- w ' ) <7 1 ) - 9 T ( a 2 + ^ 



^ ^ffs + CTiO-2 + ^ <W = 0, 

/ 2 2 4 

-4Dd e {dim e- <e (e- <e ',a 3 + dd + 



cr 2 



Here 



4)lm e -(e-,. 2 4)} 
+<r 1 Im e^V'Vi)} 



-t? T ( (7 3 + CTiCr 2 + y 



(41) 



-n 2 [Im e- 4e (e 4e ',a 2 + + a x \m e^e 1 *', <j x )'\ 



(42) 



(43) 



and we have used the abbreviations en = A ^£'^-n h cc 

and 



(a(e,n),/3(0,n)) 



2tt 



(a(0,n),/3(0,n))' = — 

Z7T 



J—oo 



2tt r + oo 



J-oo 



a(0, fi)/3(0, fi) 

g(n)d6dn, (44) 
a(fl, fi)/3(0, fi) 
^(fijdfldfi, (45) 



where 



9n a («) = 2 + fi o) - *'(n - n )] , (46) 

(n = n 0c = d). 



The ansatz (|39|) has not yet been inserted in ( |40| ) and 
(^2|) in order not to complicate further these equations. 
We shall keep in mind that these equations will have to 
be modified later when the solution of (^) is inserted 
in ( p| ) and (|4^). The linear equations in the preceding 
hierarchy should have solutions periodic in 6. A solution 
of the homogeneous equation CU — could be added to 
the solutions of the linear nonhomogeneous equations in 
the preceding hierarchy. However, all such terms are to 
be omitted, for the amplitude A(T; e) already takes care 
of them. 

The solution of @, 



A T e l6 



D + itt D + iVL 



(47) 



is 



<T2 + T = 



A 2 e 2W 



Are 1 



(D + iCL){2D + {D + ity 2 



(48) 



Note that one term proportional to e 10 /(D + in) and 
other terms decaying on the fast scale t are solutions of 
CU = and could have been added to (Q). According 
to what was said above, all such terms are to be omitted. 

We now insert fl4g| ) in ( p4| ) and (|42] ) together with the 
CE ansatz (p9). The solution of (|4ip is 



cr 3 + CTi(T 2 + — 

6 



K 2 - 4» 2 
4D(D + iST) 
A\A\ 2 



A 



(D + ifl) 2 (2D + 
AA- 1 1 L 1 



{D + m) 3 

9 +CC 



T \D+m ' 2D+m 
(D + ifl)(2D + m) 



„2i8 



+ 



A 3 



;mi> 



(D + ifi) (2D + in) (3D + in) 



cc 



cc. 



(49) 



From this we obtain ct 3 , and finally, from (f42|), cr 4 . To 
obtain the leading order approximation, wc only need 
to determine A(T; e). Now, ( J49| ) holds provided that the 
nonresonance condition (needed to remove secular terms) 



x D + itt 



—,P(Sl,T;A)) = 



(50) 



holds, where P(n, T; A) denotes the coefficient of e %B on 
the right-hand side of Q) || . Equation (^) yields 

= ^(A^ 2 ~4n 2 )A+^| 2 A = 0. (51) 
I 5 

The function is determined from a similar nonres- 
onance condition for (E3): the coefficient Q(£l,T;A) of 



G 



e on the right-hand side of ([42) should also satisfy (50) 



The result is 



Kn 



'-An 



(\A\ 2 A) : 
5D 



23 
25D 



\A\ 2 A 7 



Insertion of Equations (|5l|) and (|)2|) into (| 
sought amplitude equation: 



(52) 



yields the 
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y(if 2 -4r> 2 )A-^| 2 A = 



i 



0(e 2 ). (53) 
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This equation is the scaled "normal form" studied by 
Dangelmayr and Knobloch in |Q [cf. their equations 
(3.3), p. 2480]. The general analysis developed in that 
Reference for general scaled normal forms was employed 
in (|] to study (|53|) and will not be repeated here. 



IV. CONCLUSIONS 



We have applied a modified Chapman-Enskog method 
to two problems related to Kuramoto models of synchro- 
nization of globally coupled phase oscillators. First of 
all, we found a consistent two-term Smoluchowski ap- 
proximate equation for a model of oscillators with iner- 
tia in the limit of small inertia (as mD — > 0+). Second, 
we modified the Chapman-Enskog method to find di- 
rectly the scaled normal form corresponding to the 0(2)- 
symmetric Takens-Bogdanov bifurcation at the tricritical 
point of a standard Kuramoto model with bimodal distri- 
bution of oscillator natural frequencies. Key ingredients 
of the CEM are: (i) solving a zeroth-order problem whose 
solution is determined up to certain amplitude functions; 
(ii) assuming an expansion for the solution all whose 
higher order terms depend on (slow) time through the 
amplitude functions only; (iii) assuming that the right 
sides of the equations of motion for the amplitude func- 
tions are expansions whose coefficients are functionals of 
the amplitudes. These coefficients are determined by ap- 
propriate solvability conditions for the hierarchy of linear 
equations resulting from insertion of all these assump- 
tions in the original equations. Collecting the desired 
number of coefficients, we obtain approximate equations 
of motion for the amplitude functions as the result of the 
method. I believe that the techniques explained in the 
present paper will be useful in many other problems of 
physical interest. 
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